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Abstract 

Recently, several research efforts showed that the analysis of joint 
spectral characteristics of sets of matrices is greatly eased when these 
matrices share an invariant cone. In this short note we prove two new 
results in this direction. 

We prove that the joint spectral subradius is continuous in the 
neighborhood of sets of matrices that leave an embedded pair of cones 
invariant. 

We show that the (averaged) maximal spectral radius, as well as 
the maximal trace, of products of length t, converge towards the joint 
spectral radius when the matrices share an invariant cone, and addi- 
tionally one of them is primitive. 

1 Introduction 

The Perron- Frobenius theorem is one of the most well known theorems in lin- 
ear algebra. It states strong properties that matrices with nonnegative entries 
enjoy. Starting in the fifties, is has become clear that not only this theorem, 
but many other properties of nonnegative matrices can be generalized to 
much more general cases, the fundamental feature being of importance for a 
matrix A is that there exists an invariant proper cone K, that is, a proper 
cone such that AK C K. (A cone K is said proper if it is closed, convex, 
with nonempty interior, and contains no straight line. Unless explicitly men- 
tioned, all cones are supposed to be proper in the following.) If K is an 
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invariant cone for A, we say that A is K-nonnegative. If AK C intfT, we say 
that A is K-positive. Finally if there exists a natural number t such that A 1 
is .fT-positive, we say that A is primitive. (All these definitions are obvious 
generalizations of the case of nonnegative matrices with K = R™.) As an 
example of the generalizations one can obtain, let us mention the generalized 
Perron-Frobenius theorem (see for instance [TU]). 

Theorem 1. Generalized Perron-Frobenius Theorem. Let K be a proper 
cone. If a matrix P is K -primitive, then it has a single eigenvalue of largest 
modulus, which moreover is a real positive number. 

The interested reader can find a survey of properties of K- nonnegative 
matrices in [18]. 

Our goal in this short note is to exploit the assumption of .fT-nonnegativity 
for the study of finitely generated matrix semigroups, and more precisely of 
joint spectral characteristics. The joint spectral characteristics of a set of 
matrices are quantities that allow to describe the asymptotic behaviour of the 
semigroup generated by this set, when the length of the products increases. 
In this paper we will restrict our attention to two of these quantities: 

Definition 1. For a set of matrices X C R nxn , the joint spectral radius p(S) 
and joint spectral subradius are respectively defined as: 

p(S) 4 li mm ax{||A 1 ---A t ir A : AeS}, (1) 

t— »oo 

p(E) = \immm{\\A ll ---A it \\ 1 ^ : A e £}. 

t— >oo 

Both these limits exist and do not depend on the norm chosen. They are 
natural generalizations of the notion of spectral radius of a matrix (i.e., the 
maximal modulus of the eigenvalues) to a set of matrices. The joint spectral 
radius appeared in [T7], and the joint spectral subradius in [8]. Not only these 
quantities account for stability issues of switching linear systems, but they 
have also found many applications in various different fields of Engineering, 
Mathematics, and Computer Science. See [9] for a recent survey on these 
quantities. 

The issue of computing the joint spectral radius has also been largely studied, 
and several negative results are available in the literature. For instance, the 
issue of exactly computing the joint spectral radius, as well as the subradius, 
are both known to be undecidable. This fact has led to a rich literature, where 
techniques from very different fields (from Control Theory to Ergodic Theory, 
Automata Theory, etc.) have been applied to the topic (see [T| \5\\7\ IT3HTB] as 
a few examples). 

The joint spectral subradius has been much less studied until recently. In 
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fact, the only methods we know of that allow to compute the joint spectral 
subradius have been proposed only very recently [BICES]. These methods are 
proved to be efficient only in some favorable cases, namely, when the set of 
matrices share a common invariant embedded pair of cones (see below for 
definitions). It appears that this assumption allows for more appealing prop- 
erties, and in this note we prove (Section [2]) that under the same hypothesis, 
the joint spectral subradius is a continuous function. Then, in Section El we 
study the joint spectral radius, and we generalize a recent result on nonneg- 
ative matrices to any set leaving a cone invariant. This result is about the 
convergence of the maximum trace, and the maximum spectral radius of the 
products towards the joint spectral radius, when the length of the products 
increases. 

The joint spectral characteristics have attracted much attention in recent 
years, especially in the case of nonnegative matrices because many appli- 
cations involve such matrices [5J [101 H21 EO] • These applications make the 
theoretical study of (K-) nonnegative matrices of particular importance. 

2 Continuity of the subradius 

The continuity of the joint spectral radius (w.r.t. the Hausdorff distance^) is 
well-known: 

Proposition 1. I21\j For any bounded set of matrices X G M n , and for any 

e > 0, there is a 5 > such that 



Less is known on the subradius. This quantity is not continuous, as shown 
on the next example, drawn from [9]: 

Example 1. The sequence of sets 



The Hausdorff distance measures the distance between sets of points in a metric space: 



D(E,E')<5^|p(E)-p(E')|<e. 




converges towards 




(£,£') = max { sup { ini l\\A-A'\\}, sup { ini ML4 - A'\\}} 
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when k — > oo. For any k G N, we have p(Efc) = because the product (AiA^) 2 
is the zero matrix. However p(E) = 1 (This is because the lower right entry 
of any product is equal to one) and thus, the joint spectral subradius is not 
continuous in the neighborhood of E. 

Thus, the nonnegativity of matrices is not sufficient for ensuring the con- 
tinuity of the subradius in the neighborhood of the set. However, we show in 
this section that if the matrices share a second invariant cone, then we have 
the continuity. More precisely, let us consider a cone K G MJ 1 . We say that 
a convex closed cone K' is embedded in K if (K' \ {0}) C inti^. In this case, 
following [TB], we call {K,K'} an embedded pair. Note that the embedded 
cone K' may be degenerate, i.e., may have an empty interior. An embedded 
pair {K,K'} is called an invariant pair for a matrix A (or a set of matrices 
E) if the cones K and K' are both invariant for A (for the matrices in E). 
The following definition aims at characterizing the "embeddedness" of the 
pair (K, K') . 

Definition 2. [TSf For a given embedded pair {K, K'} the value j3(K, K') is 
the smallest number such that for any line intersecting K and K' by segments 
[x,y] and [x',y'} respectively (with [x,x'] C [x, y'\) one has 1 < J^~^,j < (5. 

In the following we denote by E* the set of products of length t of matrices 
in E. In the developments below we use the two following results: 

Lemma 1. f77]/ Let E be a compact set of matrices that share an invariant 
cone K. If there exists a nonzero vector x G K such that 

WA G E, Ax >k rx, 

then p(S) > r. 

Theorem 2. f7g| Theorem 2.12] For any compact set E with an invariant 
pair {K, K'}, there exists a nonzero vector x G K' such that 

VA G E, Ax > K p(E)x//3, 

where (3 = f3({K, K'}). 

We are now in position to prove the main theorem of this section: 

Theorem 3. Let E be a compact set of matrices in M. nxn , and let E& be a 
sequence of sets in M. nxn that converges to E in the Hausdorff metric. If E 
leaves an embedded pair of cones invariant, then, 

p(Efc) —> p(S) as k —7- oo. 
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Proof. Let us consider a set E that leaves an embedded pair of cones (K, K') 
invariant, such that /3(K, K') = (3. If p(S) = then the joint spectral sub- 
radius is continuous at E, being upper semi continuous and nonnegative in 
general [9]. We can then suppose that p(E) > and, by scaling the set 
of matrices, we suppose that p(E) = 2 (the joint spectral subradius is an 
homogeneous function of the entries of the matrices) . 

Fix e > 0. By upper semicontinuity, we know that there exists a 5 > 
such that 

£>(£',£)<<* => p(E') < p(E) +e. 
We still have to show that for some 5' > 0, if £)(£',£) < 5', then 

p(E') > p(E) - e. 
It follows directly from the Definition fl2]) that 

p(E')=p'(E). 

Combining this with Theorem [21 there must exist an integer t and a vector 
x G K', \x\ = 1, such that 

WA e £*, Ax (2 - e/2)*x. 

Take also t large enough such that (2 - e/2)* > (2 - e)* + e/2. Since x e K' C 
int(i^), we can define rj > such that for any vector y : 

\y\<rj, =4> (e/2)x + yeK. 

Then, S being compact, there exists a 5' > such that 

£>(£',£)< 5' D(E rt , E*) < 77. 

Take such a set of matrices E', so that -D(S", E*) < 77. Then, for any product 
A' G E'*, taking A the corresponding product in E* close to A', 

A'x = {A' - A + A)x 

> K (2 - e/2fx + (A' - A)x, 

> K (2 - efx + {A' - A)x + (e/2)x 

>^ (2-e)*x. 

By Lemma dj this implies that p(E') > 2 — e. □ 

Corollary 1. The joint spectral subradius is continuous in the neighborhood 
of any set of matrices with positive entries. 

Proof. It is known (see [161 Corollary 2.14]) that if all matrices of a set E 
are positive and in each column of any matrix the ratio between the greatest 
and the smallest elements does not exceed c, then E has an invariant pair of 
embedded cones, for which (3 < c 2 . □ 
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3 Asymptotic regularity 



It is well known that for any bounded set of matrices, the joint spectral radius 
can be alternatively defined in terms of the maximal spectral radius (instead 
of the maximal norm): 

Theorem 4. (Joint Spectral Radius Theorem, f^j) 
For any bounded set of matrices E, 

limsupsupV^'Oi) : A E E*} = p(E). 

t— >oo 

Some attention has been recently given to a peculiarity in the alternative 
definition provided by Theorem HJ the maximal averaged spectral radius 
asymptotically converges towards the value of the joint spectral radius in 
limit superior, but not at every time steps. Simple examples illustrate this 
fact, as for instance 

o or ll 0. 



(2) 



For this set of matrices p(S) = 1 but all products of odd length have a 
spectral radius equal to zero. It is thus natural to try to understand when 
the maximal averaged spectral radius actually converges or to give sufficient 
conditions for it. 

Moreover, the next proposition shows that the joint spectral radius can 
also be defined as the limit superior of the rate of growth of the maximal 
trace in the semigroup: 

Proposition 2. |7|/ For any finite set of matrices, the joint spectral radius 
satisfies 

p(S) = limsupmax {trace^iA)}. (3) 

Here again, the joint spectral radius is defined as a limit superior, and one 
could wonder when the sequence actually converges. A sufficient condition 
relying on the nonnegativity of the matrices has been proposed recently: 

Theorem 5. l2Ef Let E = {A\, . . . , A n } be a set of nonnegative matrices. If 
one of them is primitive, then 

max {trace{ A) 1,1 } ->• p(E) (4) 

as t — >■ oo. 

In this section we generalize Theorem to arbitrary invariant cones: 
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Theorem 6. Let a set of matrices X = {Ai,...,A m } share an invariant 
cone K. If there is a matrix £ S which is i^-primitive, then, both the 
quantity 

max {trace l/t { A) j 
AeE* 

and 

m&x{p 1/t (A)} 
converge towards p when t tends to oo. 

The proof is inspired from developments in [22] and results from (3] which 
we now recall. 

Definition 3. Kronecker product. Let A £ R niXn \ B £ W^ xm \ The 

Kronecker product of A and B is a matrix in pi" 1 ^^^ defined as 

/ A hl B ... A hn2 B\ 
(.1, B) A : : : 

\A ni iB ... A ni ^ n2 B J 

The k-th Kronecker power of A, denoted A® k , is defined inductively as 

A® k = A^A®^^ A m = A. 
Proposition 3. 1. For any matrices A,B,C,D, 

{A ®B)(C®D) = {AC) <g> (BD). 

2. For any matrix A, we have 

trace(A m ) = trace k (A). 

The main reason for introducing Kronecker products is the following the- 
orem: 

Theorem 7. |2J/ Let X = {Ax, . . . ,A m } be a set of matrices that leaves a 
cone invariant. Then, we have the following property: 

p(E) = lim p 1 ' k (Af k + ■■■ + A® k ). (5) 

fc— >oo 

Before to prove Theorem we prove a last technical lemma that will be 
necessary in the proof. 
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Lemma 2. Let K be a proper cone. If a matrix A is K -primitive, then for 
any k, A® k is K® k -primitive, where 

K ®k a Conv{xi <g • • • <g x k : x t G K, i = l,...,k} 

is a proper cone. 

Proof. The fact that K® k is a proper cone is proved in [31 Lemma 4]. 

It is clear that if A is .fT-nonnegative, then A® k is IT® -nonnegative, since 

A m (x 1 <g • - • <g x k ) = Axi ® ■ ■ • <g Ax k . 

It remains to show that all extremal points of K® k are mapped in the interior 
of K® k . 

So, let us take an extremal point x = X\ <g • • • <g x k of K® k , and show that 
for some t > 0, (A® fc )'x G intJC® fc . Since 

(A® fc )'x = A t x 1 <g • • • <g A*x fc 

and A is primitive, it suffices to show that 

Xi G intK => Xi <g • • • <g x fc G intftT® fc . 

We will in fact show the slightly stronger property that for any set of proper 
cones Ki, 

Xi G inti^i =>- Xi <g • • • <g Xk G int(Ki <g • • • <g 

where 

K i <g) • • • (g) = Conv{s! (g ■ ■ ■ (g a;^ : Xj G i^}. 

By induction, it is sufficient to prove it for k = 2. Let us be given two proper 
cones K\ G R™ 1 , K 2 G IR™ 2 , Xi G int(ifj). Suppose by contradiction that 
Xi<gx 2 G ^(gfGVnt^itgiG)- Then, there exists a vector z G R™ lX?12 , z ^ 0, 
such that (z,xi (g x 2 ) = 0, (z,w) > OWw G If i <g if 2 . Now, since for any i, 
1 < i < rii, there exists a A small enough such that (xi + Ae^) <gx 2 G K 1 ^)K 2 
(e, is the ith standard basis vector); this implies that 

(z, (xi + Aej)) (g x 2 = (z, xi (g x 2 ) + (z, Xei <g x 2 ) = (z, Ae^ (g x 2 ) > 0, 

and thus (z, ej (g x 2 ) = (because A can be positive and negative). For the 
same reason, for any j, (z, x\ (g ej) = 0. Now, 

(z, (xi + \ei) (g (x 2 + 7ej) = {z, xi <g x 2 ) + (z, Ae^ (g x 2 ) 

+(z, xi <g jej) + (z, \^Ci <g ej) 
= (z, \^Ci <g ej) 
> 0. 

Thus, (z, ei®ej) = for all i,j, which implies that z = 0, a contradiction. □ 
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We are now in position to prove the main result of this section. 

Proof. (Proof of Theorem [6]) We first claim that for any natural number k, 
the kth Kronecker powers of the matrices in E satisfy the property 

lim trace^Af + • ■ ■ + A^) 1 = p(Af k + ■■■ + A® k ). 

t—>oo 



By Lemma [2] above, one of the matrices Af k is i^® fc -primitive. Moreover, it 
is well known that a sum of nonnegative matrices, one of which is primitive, 
is actually also primitive, thus, so is 



Finally, by Theorem [TJ it is clear that 



lim(trace(P*) lA ) = p{P) 



t— >oo 



and the claim is proved. 

Now, remark that for any k, t 6 N, 



trace((^ A® k f) = ^ trace(A® fc ) < m'(maxtrace fc (A)). (6) 



AGE' 

Age age* 



For deriving the above relations, we successively used items 1 and 2 of Propo- 
sition [3j Dividing this equation by m* and taking to the power l/(kt), we 
obtain 

trace 1 /( tk )((22A 9h /m) t ) < max {trace 1 /* (A)}. (7) 
Age AeS 

Now, for any fceN, the lefthand side in (JZJ) tends towards P (Sags A® k /m) 
as t — > oo. 

This proves part one of the theorem, as the latter quantity is arbitrarily close 
to for large k (Lemma [7]). 

For the second part of the theorem, since p(A) > tr&,ceA/n (n is the 
dimension of the matrices), we have 



max {(trace (A) /n) 1,1 } < max {p(A) 1 }, 

Age* Age' 

which immediately implies that max J 4 6 E'{/0 1 '*C<4-)} — > p when t — > oo. 



□ 
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